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1.  INTRODUCTION 


An  important  problem  with  adaptive  arrays  [1,2]  is  that  their  performance  de¬ 
teriorates  with  interference  bandwidth.  The  wider  the  bandwidth  of  an  interference 
signal,  the  more  difficult  it  is  for  an  adaptive  array  to  null  it  [3, 4, 5, 6, 7). 

Figure  1.1,  from  [7],  illustrates  this  problem.  It  shows  the  output  SINR  (desired 
signal-to-interference-plus-noise  ratio)  achieved  by  a  two  element  adaptive  array 
when  an  interference  signal  arrives  from  angle  0,.  The  figure  assumes  two  isotropic 
elements  a  half-wavelength  apart,  a  desired  signal  with  a  0  dB  SNR  (signal-to-noise 
ratio  per  element)  arriving  from  broadside  (0j  =  0°),  and  interference  with  a  40  dB 
INR  (interference- to-noise  ratio  per  element).  The  SINR  is  shown  for  several  values 
of  B,  where  B  is  the  relative  bandwidth,  i.e.,  the  ratio  of  the  absolute  bandwidth  to 
the  center  frequency1.  As  may  be  seen,  for  B  =  0.02  the  output  SINR  has  dropped 
about  3  dB  below  its  value  for  B  =  0.  Larger  bandwidths  cause  increasingly  more 
degradation. 

When  the  bandwidth  performance  of  an  adaptive  array  is  inadequate,  three 
methods  exist  for  improving  its  performance: 

1.  Using  more  elements  in  the  array, 

2.  Using  tapped  delay-lines  behind  the  elements,  and 

3.  Using  FFTs  (Fast  Fourier  Transforms)  behind  the  elements. 

Let  us  compare  the  performance  obtained  with  each  of  these  approaches. 


First,  suppose  we  add  more  elements  to  the  array.  Figure  1.2  shows  the  SINR 


when  the  array  has  3,  5,  10  or  20  elements,  instead  of  2  as  in  Figure  1.1.  (In 
each  case  the  array  elements  lie  along  a  straight  line  with  a  half  wavelength  spacing 
between  elements.  All  other  parameters  are  the  same  as  in  Figure  1.1.)  For  each 
array  size,  the  SINR  is  shown  for  B  =  0  and  B  =  0.2.  As  may  be  seen  from  Figures 
1.1  and  1.2,  adding  elements  does  improve  the  SINR.  However,  it  is  interesting  that 
no  matter  how  many  elements  are  used,  there  is  always  a  region  for  9,  near  where 
the  SINR  for  B  =  0.2  is  poorer  than  for  5=0. 

The  second  way  to  improve  bandwidth  performance  is  to  use  a  tapped  delay-line 
behind  each  element.  (Figures  1.1  and  1.2  assume  a  single  complex  weight  behind 
each  element.)  Figure  1.3  shows  a  two-element  array  with  a  2-tap  delay  line  behind 
each  element.  Figure  1.4  shows  the  output  SINR  vs.  0,  achieved  by  this  array  when 
B  =  0.2  and  when  each  delay  is  a  quarter  wavelength  long  at  the  carrier  frequency. 
(The  other  parameters  are  the  same  as  in  Figure  1.)  Comparing  Figure  1.4  with 
Figure  1.1  shows  that  the  delay  lines  behind  the  elements  have  fully  overcome  the 
bandwidth  degradation.  The  performance  for  B  =  0.2  in  Figure  1.4  is  just  as  good 
as  that  for  B  =  0  in  Figure  1.1. 

The  third  method  that  has  been  proposed  for  improving  bandwidth  performance 
is  to  use  an  FFT  (Fast  Fourier  Transform)  behind  each  element  of  the  array,  with 
a  separate  weight  on  each  frequency  bin.  This  approach  is  shown  in  Figure  1.5.  At 
the  output  of  each  element,  an  A/D  converter  takes  samples  of  the  received  signal. 
When  K  samples  are  available  from  each  element,  these  samples  are  transformed 
with  an  FFT.  The  FFT  produces  K  frequency  domain  samples  of  each  element 
output.  These  frequency  domain  samples  are  each  multiplied  by  a  weight  and  then 
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added  to  the  corresponding  frequency  bins  from  other  elements.  Finally,  the  inverse 
FFT  (IFFT)  is  taken  of  the  frequency  domain  samples  to  obtain  the  time-domain 
array  output  samples. 

The  use  of  FFTs  behind  the  elements  has  a  certain  intuitive  appeal  as  a  method 
of  improving  array  bandwidth  performance.  In  effect,  the  FFT  divides  the  signal 
bandwidth  into  smaller  subbands.  (For  this  reason,  this  technique  is  sometimes 
called  frequency  subbanding.)  With  a  separate  weight  on  each  frequency  bin,  the 
array  can  compensate  differently  for  each  frequency  subband. 

Unfortunately,  calculations  of  array  output  SINR  for  FFT  processing  often  in¬ 
dicate  poorer  performance  than  that  for  tapped  delay  lines.  Figure  1.6  shows  a 
typical  set  of  results.  It  shows  the  output  SINR  from  the  same  two-element  array 
as  in  Figure  1.1,  but  with  a  K-point  FFT  behind  each  element.  The  curves  are 
computed  for  the  same  bandwidth  of  B  =  0.2  and  for  a  sampling  interval  that 
makes  one  period  of  the  FFT  frequency  response  equal  to  the  signal  bandwidth. 
The  SINR  is  shown  for  K  =  2,4,8  and  16  samples  in  the  FFTs.  As  may  be  seen, 
the  performance  does  improve  with  K ,  but  even  for  K  =  16  it  is  not  as  good  as  the 
performance  for  tapped  delay  lines  with  only  two  taps,  as  seen  in  Figure  1.3.  Thus, 
in  spite  of  the  intuitive  appeal  of  FFTs,  their  performance  can  be  disappointing. 

The  present  study  was  done  in  an  effort  to  understand  the  relationship  between 
tapped  delay-line  and  FFT  processing  in  adaptive  arrays.  The  study  was  motivated 
by  the  fact  that  there  seemed  to  be  no  apparent  reason  why  FFT  performance  should 
be  poorer  than  tapped  delay-line  performance. 

In  this  report  we  show  that  inserting  an  FFT  between  the  delay-line  taps  and 
the  weights  in  an  adaptive  array  in  fact  has  no  effect  on  the  output  SINR.  The 
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difference  in  performance  between  the  two  approaches  noted  above  is  simply  due 
to  the  different  time  delays  between  samples  or  different  numbers  of  samples  used 
in  each  case.  When  a  tapped  delay-line  array  and  an  FFT  array  use  the  same  time 
between  samples  and  the  same  number  of  samples,  their  performance  is  identical. 

We  proceed  as  follows.  In  Section  2,  we  first  prove  that  inserting  any  linear, 
invertible  transformation  between  the  delay-line  taps  and  the  weights  in  an  adaptive 
array  has  no  effect  on  the  array  output  SiNR.  We  prove  this  theorem  for  both 
the  LMS  array  and  the  Applebaum  array.  (For  the  Applebaum  array,  the  proof 
holds  only  if  the  steering  vector  is  also  transformed  in  the  appropriate  way.)  Next, 
in  Section  3,  we  show  that  FFTs  in  an  array  simply  provide  a  linear  invertible 
transformation  of  the  signals.  Taken  together,  these  two  results  show  that  FFT 
processing  is  entirely  equivalent  to  tapped  delay-line  processing.  Then,  in  Section  4, 
we  discuss  how  FFT  processing  affects  the  array  weights  and  the  covariance  matrix 
eigenvalue  spread.  We  also  discuss  how  tapped  delay-line  and  FFT  parameters 
are  often  chosen,  and  point  out  how  these  choices  result  in  performance  differences 
between  the  two  methods.  Finally,  Section  5  contains  the  conclusions. 
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2.  A  SIMPLE  THEOREM 

In  this  section  we  present  a  simple  theorem  about  adaptive  arrays  with  tapped 
delay  line  processing.  We  show  that  inserting  an  arbitrary  invertible  linear  trans¬ 
formation  between  the  delay  line  taps  and  the  adaptive  weights  has  no  effect  on 
either  the  output  signal  or  the  output  SINR.  We  prove  this  result  for  LMS  arrays  in 
Part  2.1  and  for  Applebaum  arrays  in  Part  2.2.  For  Applebaum  arrays,  the  steering 
vector  must  be  transformed  in  the  proper  way  along  with  the  signals  for  the  theo¬ 
rem  to  hold.  We  use  this  theorem  in  the  next  section  to  show  that  FFT  processing 
is  equivalent  to  tapped  delay-line  processing. 

2.1  The  LMS  Array 

Consider  an  adaptive  array  with  M  elements,  as  shown  in  Figure  2.1.  Assume 
each  element  is  followed  by  a  tapped  delay-line  with  K  taps  and  a  delay  of  To  seconds 
between  taps.  The  output  of  the  first  tap  behind  each  element  is  the  element  signal 
itself,  with  no  delay.  Let  xm*(£)  denote  the  (analytic)  signal  from  element  m  at  tap 
k.  Then  xmi(£)  is  the  signal  received  on  element  m,  and 


imk(t )  =  xml(£  -  \k  -  1]T0). 


(2.1) 


Let  us  suppose  the  xm*(<)  are  used  as  inputs  to  an  LMS  adaptive  array  processor 
[lj.  This  processor  multiplies  each  xm*(<)  by  a  complex  weight  wmk  and  then  adds 
the  weighted  signals  to  produce  the  array  output  sx(£),  as  shown  in  Figure  2.1. 
The  processor  uses  LMS  feedback  loops  [lj  or  an  equivalent  technique  [8]  to  set 
the  weights  to  their  optimal  values.  The  LMS  optimal  weights  minimize  the  mean- 
square  difference  between  the  array  output  .s T(t)  and  a  reference  signal  r(t),  as 
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Figure  2.1:  An  AZ-element  Array 

described  in  [lj.  Moreover,  as  long  as  the  reference  signal  is  correlated  with  the 
desired  signal  and  uncorrelated  with  interference,  the  optimal  weights  also  maximize 
SINR  at  the  array  output  [9] . 

For  a  given  set  of  tap  signals  xmk(t),  the  optimal  weights  may  be  calculated  as 
follows.  Let  Xm(t)  and  Wm  (with  1  <  m  <  M)  be  column  vectors  containing  the 
signals  and  weights  at  the  K  taps  behind  element  m,  i.e., 

Xm{t)  =  xm2(t),  ...,  xmK{t)}T,  (2-2) 

and 

wm  =  [Wml,  Wm 2,  WmK\T .  (2.3) 

(Superscript  T  denotes  transpose.)  We  refer  to  .Ym(t)  as  the  element  signal  vector 
and  to  Wm  as  the  element  weight  vector.  Then  let.  X(t)  and  U  be  the  total  signal 
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st  =  E\X"r{t)  j. 


(2.8) 


In  these  equations,  *  denotes  complex  conjugate  and  r(t)  is  the  reference  signal 
[l].  We  assume  the  signals  xmk{t)  as  well  as  r(t)  are  all  jointly  stationary  random 
processes,  so  $x  and  Sx  do  not  depend  on  t. 

With  the  weight  vector  W  given  by  (2.6),  the  array  output  signal  is 

3,(0  =  XT{t)W  =  XT{t)<t>-xlSx.  (2.9) 


Now  consider  the  following  alternative  situation.  Suppose  that,  instead  of  using 
the  signals  x mk{t)  as  inputs  to  the  processor,  we  first  combine  them  in  some  manner 
to  produce  a  new  set  of  signals  ymk{t),  where  1  <  m  <  M  and  1  <  k  <  K. 
Specifically,  suppose  Ym  is  the  /^-component  vector 


Xm  —  (yml(O)  ym2(0i  •••5  i/m/f(0]  > 


and  Y  is  the  M/f-component  vector, 

f  n  1 


Y(t)  = 


Ym 


Then  let  us  assume  that 


(2.10) 


(2.11) 


rm  =  txu). 


(2.12) 


TRANSFORMATION  T 


WEIGHT 

CONTROL 

SYSTEM 


ARRAY  / 
PROCESSOR 


Figure  2.2:  Tapped  Delay-Line  Array  with  Transformation  T 

where  T  is  an  MK  x  MK  matrix  of  constants.  Thus,  each  ymk{t)  is  a  linear 
combination  of  the  xmk(t).  Now  let  us  use  the  ymk(t)  as  inputs  to  the  same  LMS 
processor  as  before.  Figure  2.2  shows  the  new  arrangement  with  transformation  T 
between  the  xmk[t)  and  the  ymk(t).  For  this  case,  we  denote  the  mkth  array  weight 
by  umk,  to  distinguish  it  from  wmk  in  Figure  2.1. 

With  the  signals  ymk{t)  as  inputs,  the  LMS  processor  will  produce  optimal 


weights  given  by 


U  = 


where  U  is  the  new  weight  vector, 


U  —  [u  1 X ,  U 12 1  •  •  •  1  U21i  •  •  •  1  “mk! 


(2.13) 


(2.14) 


is  the  covariance  matrix  associated  with  the  signals  ymk{t), 


*y  =  E[Y'(t)YT(t)\, 


and  Sv  is  the  reference  correlation  vector  for  the  signals 


(2.15) 


5,  =  £[y*(*)r(t)]. 


(2.16) 


f(t)  is  the  same  reference  signal  as  in  (2.8).  The  array  output  signal  for  this  case, 


sv(t),  is 


3y(t)  =  YT(t)[/  =  yT(t)*rls¥. 


(2.17) 


Now  it  is  easy  to  show  that  if  T  is  invertible,  s„(t)  and  sz(t)  are  in  fact  identical 
signals.  Substituting  (2.12)  into  (2.15)  and  (2.16),  and  using  (2.7)  and  (2.8),  we 


<f>v  =  E[Y-{t)YT{t)\  =  E\T'  X'  (t)  XT  (t)TT\  =  T’$zTt, 


(2.18) 


Sv  =  E[Y'{t)r{t)\  =  E[T'X'(t)f{t)\  =  T'SX. 

If  T  is  invertible  (i.e.,  nonsingular),  s„(t)  in  (2.17)  reduces  to 

mo  = 

=  XT{t)TT[T'$zTT}-lrSz 
=  xT(t)TT[TT}-l<t>;1[T‘}-lrsz 
=  xT(t)<trzlsz 

=  MO- 


(2.19) 


(2.20) 


Furthermore,  the  output  SINR  is  the  same  for  the  two  arrays.  Substituting 
(2.18)  and  (2.19)  into  (2.13)  shows  that  the  weight  vectors  U  and  W  are  related  by 


U  =  \T^xTt]~1TSx 
=  [Tt]-1*-1[T'\-1T'Ss 
=  \TT)-l^Sz 

=  [Tt}-1W.  (2.21) 

Consider,  for  example,  the  output  desired  signal.  Suppose  Xd(t)  is  the  desired  signal 
part  of  signal  vector  X(t)  and  Yd(t)  is  the  desired  signal  part  of  F(t).  Then 


Y&)  =  TXd(t).  (2.22) 

The  array  output  desired  signal  in  Figure  2.1  is 

s.,(0  =Xj(t)W,  (2.23) 

whereas  the  output  desired  signal  in  Figure  2.2  is 

*«At)  =  YlU.  (2.24) 

However,  using  (2.21)  and  (2.22)  in  (2.24)  gives 

3W(0  =  [TXd(t))T[TT\-lW  =  (2.25) 


Thus,  the  arrays  in  Figures  2.1  and  2.2  have  identical  output  desired  signals,  and 
hence  identical  output  desired  signal  powers.  A  similar  argument  shows  that  the 
output  interference  and  thermal  noise  powers  are  also  identical  for  the  two  arrays. 
From  this  it  follows  that  the  output  SINR  is  the  same  in  Figures  2.1  and  2.2. 


Thus,  inserting  a  linear  transformation  T  between  the  elements  and  the  proces¬ 
sor,  as  in  Figure  2.2,  has  no  effect  on  the  output  signal  or  the  output  SINR.  The 
only  requirement  for  this  result  to  hold  is  that  the  transformation  be  invertible. 

2.2  The  Applebaum  Array 

Now  suppose  the  array  processor  in  Figure  2.1  uses  Applebaum  control  loops  [2] 
or  an  equivalent  technique  [8]  to  set  the  weights  wmk-  Applebaum  control  loops  use 
a  steering  vector  instead  of  a  reference  signal  to  maintain  the  array  response  in  the 
desired  signal  direction.  With  Applebaum  loops,  the  steady-state  weight  vector  in 
the  array  will  be  [2] 

W  =  (2.26) 

where  /i  is  an  arbitrary  (nonzero)  gain  constant  and  V  is  the  steering  vector.  The 
array  output  signal  will  be 

3,(0  =  XTW  =  nXr{t)*-xlV.  (2.27) 

Now  suppose  a  transformation  T  is  inserted  between  the  signals  xmk(t)  and  the 
Applebaum  processor,  as  in  Figure  2.2.  Let  the  processor  now  have  a  new  steering 
vector  Q.  The  steady-state  weight  vector  in  this  case  is 

U  =  (2.28) 

If  (2.18)  is  substituted  for  (2.28)  becomes 

U  =  n[TT\-l*:l\r\-lQ,  (2.29) 
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so  the  array  output  signal  is 

i,(t)  =  YT{t)U 

=  vXT(t)TTlTT]-l4>-xlir}-1Q 

=  nXTm;1[r)~1Q-  (2-30) 


Comparing  (2.30)  with  (2.27)  shows  that  sz[t)  and  sy(t)  will  be  the  same  if 


V  =  {T'}-lQ, 


(2.31) 


s*  V  v* 


i.e.,  if 

Q  =  TV.  (2.32) 

Thus,  inserting  the  transformation  T  between  the  elements  and  the  processor  has 
no  effect  on  the  output  signal  (or,  as  with  the  LMS  array,  on  the  output  SINR),  if 
the  steering  vector  is  transformed  according  to  (2.32). 

We  have  now  shown  the  major  result  of  this  section:  placing  an  invertible  trans¬ 
formation  T  between  the  delay-line  taps  and  the  adaptive  processor  has  no  effect 
on  the  output  signal  or  the  output  SINR.  This  result  holds  for  the  LMS  array  and 
also  for  the  Applebaum  array  if  the  steering  vector  is  properly  transformed. 

It  is  important  to  note  what  this  result  says  about  array  bandwidth  perfor¬ 
mance.  Since  the  transformation  T  has  no  effect  on  the  output  signal  or  SINR  of  an 
array,  it  also  can  have  no  effect  on  the  bandwidth  performance  of  an  array.  In  other 
words,  the  theorem  applies  no  matter  what  signals  are  present  in  the  array.  What¬ 
ever  the  signal  bandwidths  are,  the  array  SINR  will  be  the  same  with  or  without 
the  transformation  T.  Thus,  there  is  no  invertible  transformation  T  that  one  can 
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3.  AN  ARRAY  WITH  FFT  PROCESSING 


Now  consider  an  array  with  FFT  processing.  Such  an  array  was  shown  in  Figure 
1.5.  An  A/D  converter  behind  each  element  samples  the  signal  from  that  element 
every  T,  seconds.  The  samples  from  each  element  are  collected  in  the  input  buffer 
of  a  K-point  FFT  [10].  When  K  samples  have  been  stored,  the  FFT  is  taken.  This 
process  generates  K  frequency  domain  samples  from  each  element.  These  samples 
are  multiplied  by  weights  and  then  added  in  corresponding  frequency  bins  to  the 
weighted  samples  from  other  elements.  The  result  is  a  set  of  K  frequency  domain 
samples  of  the  array  output.  Finally,  the  inverse  FFT  (IFFT)  is  taken  to  obtain  K 
time  domain  samples  of  the  array  output.  This  entire  process  is  repeated  every  K 
samples. 

The  process  described  above  is  called  block  processing,  since  the  input  time 
samples  are  handled  in  blocks  of  K  samples.  After  each  block  of  K  input  samples  is 
collected,  a  block  of  K  array  output  time  samples  is  generated  by  the  inverse  FFT. 
Each  FFT  cycle  involves  a  block  of  K  entirely  new  samples. 

The  FFT  processing  can  also  be  done  in  a  sliding  window  mode.  In  this  case,  the 
FFTs  are  recomputed  after  each  new  time  sample,  using  always  the  most  recent  I\ 
time  samples  in  each  FFT  input  buffer.  As  we  will  see  below,  with  this  approach  it 
is  not  necessary  to  do  the  inverse  FFT  to  obtain  the  time  domain  array  output.  The 
time  domain  output  is  simply  the  sum  of  the  weighted  frequency  domain  samples. 

Block  processing  has  the  advantage  over  sliding  window  processing  that  the 
FFTs  need  be  computed  only  once  per  block,  instead  of  once  per  sample.  However, 
sliding  window  processing  has  the  advantage  that  no  inverse  FFT  is  required  to 
obtain  the  array  output.  In  this  section  we  shall  consider  both  forms  of  processing. 


il 


We  first  define  notation  for  the  sampled  element  signals.  To  make  the  notation 
here  consistent  with  that  in  Section  2,  we  denote  the  signal  on  the  mth  element  of 
the  array  by  xmi(t).  Let  us  concentrate  on  a  particular  set  of  K  contiguous  samples 
in  each  FFT  input  buffer  in  Figure  1.5.  Suppose  that,  of  these  K  samples,  the  most 
recent  was  taken  at  t  =  t0.  Then  the  latest  sample  of  xmi(t)  in  the  mth  FFT  input 
buffer  is  xml(t<,).  The  previous  sample  in  that  buffer  is  xmi (ta  -  Ts),  and  the  earliest 
sample  is  xmi(t0  —  [K  —  ljr,).  If  we  define  the  signal  xmk{t)  as1 

XmJc(t)  =  xmi(t  -  [k  -  l]r5),  (3.1) 

we  may  write  these  samples  as 

Xml(fo)  —  Xmi(t0), 
iml{^o  ~~  Ts )  =  Xm2(t0 ), 


- 
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Xml(fo  ~  \K  l]T»)  — 

Now  consider  the  FFT  obtained  from  these  samples.  Let  the  K  frequency  do¬ 
main  samples  produced  by  the  FFT  behind  element  m  be  denoted  by  ym i,  ym 2,  •  •  • ,  ifmK- 
These  ymn  are  given  by2 

*Eq.  (3.1)  has  the  same  form  as  Eq.  (2.1)  of  Section  2,  with  71>  replaced  by  T, . 

2  It  is  common  in  the  FFT  literature  [  10)  to  denote  the  time  domain  samples  by  Xo,  Xi,  •  •  • ,  xt<  -  i 
and  the  frequency  domain  samples  by  Xu,  Xj ,  . . . ,  X/v-  _  i .  In  this  case  the  FFT  is  usually  written 

K-L 

Xn  =  Y,  x*EKn’  o  <  n  <  K  -  1. 

fc=o 

and  the  IFFT  is 


i 

=  kE  X»EKkn’  0  <  k  <  K  —  l. 


W 


22 


■.vv.v. 


Vmn  =  £  xmk{to )4*-1)("-1),  1  <  n  <  Kt 


where 


Ek  =  e_>*. 


The  array  processor  multiplies  each  frequency  domain  sample  ymn  by  a  complex 
weight  umn.  We  assume  these  weights  are  set  to  their  optimum  values,  which  max¬ 
imize  SINR  at  the  array  output.  (The  weights  can  be  controlled  with  an  LMS  pro¬ 
cessor,  an  Applebaum  processor,  or  any  other  equivalent  processor.)  The  weighted 
samples  are  then  combined,  in  one  of  two  ways,  to  produce  the  array  output.  The 
method  used  to  combine  the  samples  depends  on  whether  block  processing  or  sliding 
window  processing  is  used. 

If  block  processing  is  used,  the  weighted  frequency  domain  samples  from  each 
element  are  added  in  corresponding  frequency  bins,  as  shown  in  Figure  1.5.  This 
step  produces  array  output  frequency  domain  samples  given  by 

M 

~fn  =  £  Umnymn.  (3.5) 

m=  1 

Then  the  inverse  FFT  of  the  /„  is  taken  to  obtain  the  time  domain  samples  of  the 
array  output.  If  we  denote  the  array  output  signal  by  s(t),  and  its  samples  by  st(f), 

h(t)  =  s(t0  -[k-  1  \T,),  1  <k<K,  (3.6) 

then  the  inverse  FFT  of  the  f„  produces  the  following  time  samples  of  s(<), 

M<o)  =  ^  £  /^K(fc“1)(n_1),  l  <  k  <  K.  (3.7) 

“  n- 1 

However,  to  make  our  FFT  notation  correspond  to  that  in  Section  2.  we  instead  write  the  FFT 
as  in  (3.3)  and  allow  the  indices  k,  n  to  vary  from  1  to  ft 
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In  a  practical  array,  the  factor  1  / K  in  (3.7}  may  be  omitted.  This  factor  is  simply  a 
gain  constant  in  the  array  output  signal,  and  it  has  no  effect  on  the  output  SINR3. 
Hence  we  assume  the  array  output  is  actually  obtained  from 

MM  =  E  fnE-K(k-1](n-l),  1  <  k  <  K.  (3.8) 

n=l 

Thus,  with  block  processing,  a  block  of  K  input  samples  is  used  in  this  manner 
to  produce  a  block  of  K  array  output  time  samples.  The  entire  process  is  repeated 
every  KT $  seconds,  using  for  each  cycle  an  entirely  new  set  of  K  samples  from  each 
element. 

If  sliding  window  processing  is  used,  on  the  other  hand,  it  is  not  necessary  to 
perform  the  IFFT  in  (3.8).  Note  that  the  most  recent  sample  of  s(t),  s(t0),  is  given 
by  (3.8)  with  k  =  1, 

K 

*(to)  =  3|(*o)  =  £/»•  (3-9) 

n—l 

Thus  s(t0)  is  just  the  sum  of  the  /„.  With  sliding  window  processing,  the  other 
samples  of  s(t)  that  could  be  found  from  (3.8)  are  not  needed,  because  an  entire 
FFT  cycle  is  performed  for  each  new  input  time  sample.  Successive  samples  of  the 
array  output  are  obtained  simply  by  repeating  (3.9)  at  each  sample  time.  Hence 
the  total  array  processing  in  this  case  is  as  shown  in  Figure  3.1. 

Note  that  the  optimal  weights  umn  in  the  processor  are  the  same  regardless  of 
whether  block  processing  or  sliding  window  processing  is  used.  The  optimal  weights 
maximize  SINR,  which  is  the  quantity, 


SINR=  p~r,  (3.io) 

where  Pd,  P,  and  Pn  are  the  average  desired,  interference  and  thermal  noise  powers 
at  the  array  output,  respectively.  If  Sd(t,),  s,(fy)  and  sn(tj)  denote  the  desired, 
interference  and  thermal  noise  signals  at  the  array  output  at  a  particular  sample 
time  tj,  then  these  powers  are  given  by 

Pd  =  £[IM<>)|2j,  (3.ii) 

P,  =  £[|M*,)I2!,  (3.12) 

and 

Pn  =  £[|3n(<,)|2].  (3.13) 

Because  we  assume  the  signals  xm*(t)  are  stationary,  the  array  output  time  sequence 
is  a  stationary  random  sequence.  Hence  each  of  the  powers  in  (3.11)-(3.13)  will  have 
the  same  value  regardless  of  which  sample  time  it  is  computed  at.  Therefore  the 
weights  that  maximize  SINR  at  one  sample  time  also  maximize  it  at  every  other 
sample  time. 

To  see  that  the  optimal  weights  with  sliding  window  processing  are  the  same 
as  those  with  block  processing,  it  suffices  to  note  that  the  array  output  for  sliding 
window  processing  in  (3.9)  is  the  same  as  the  array  output  for  block  processing  in 
(3.8)  when  k  =  1.  Hence  the  weights  that  maximize  the  SINR  at  the  k  =  1  sample 
for  block  processing  also  maximize  it  for  sliding  window  processing. 

Since  the  optimal  weights  are  the  same  for  either  type  of  processing,  and  since 
these  weights  produce  the  same  output  SINR  in  either  case,  we  shall  simplify  the 
discussion  below  by  considering  only  sliding  window  processing. 
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Figure  3.2:  An  Equivalent  Tapped  Delay-Line  Array 


Now  let  us  consider  the  relationship  between  FFT  processing  and  tapped  delay¬ 
line  processing.  First,  we  note  that  Eq.(3.l)  for  the  time  samples  in  the  FFT 
processor  has  the  same  form  as  Eq.(2.1)  for  the  signals  in  a  tapped  delay-line  pro¬ 
cessor,  except  that  the  intertap  delay  T0  in  (2.1)  is  replaced  by  the  sampling  time 
T,  in  (3.1).  Hence,  for  mathematical  purposes,  we  may  view  the  samples  in  the 
FFT  processor  as  having  been  obtained  from  tapped  delay-lines  as  shown  in  Figure 
3.2.  If  the  delay  between  taps  in  Figure  3.2  is  T and  if  every  tap  is  sampled  simul¬ 
taneously  at  t  =  to,  the  same  set  of  K  samples  will  be  obtained  from  the  tapped 
delay-lines  as  from  a  single  A/D  converter  behind  each  element. 
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Second,  we  note  that  the  frequency  domain  samples  ymn  are  each  a  linear  com¬ 
bination  of  the  input  samples  im*(to)-  The  linear  combination  is  just  the  FFT  in 
(3.3).  Third,  we  have  shown  that  the  array  output  (for  sliding  window  processing) 
is  just  the  sum  of  the  weighted  frequency  domain  samples  as  in  (3.9).  Hence,  the 
array  with  FFT  processing  is  mathematically  equivalent  to  an  array  with  tapped 
delay-lines,  followed  by  a  linear  transformation  of  the  signals,  followed  by  weighting 
and  summing,  as  shown  in  Figure  3.2. 

Moreover,  note  that  the  A/D  converters  at  the  delay-line  taps  in  Figure  3.2 
play  no  fundamental  role  in  the  operation  of  the  array.  The  same  array  output 
samples  would  be  produced  by  eliminating  the  A/D  converters  in  Figure  3.2  and 
instead  putting  a  single  A/D  converter  at  the  array  output  as  in  Figure  3.3.  The 
array  is  then  an  analog  adaptive  array  with  tapped  delay-lines,  followed  by  an  A/D 
converter  at  the  array  output.  The  A/D  converter  in  Figure  3.3  serves  only  to 
discretize  the  array  output,  but  has  no  effect  on  the  output  SINR  of  the  array. 

The  transformation  between  the  xm*(t0)  and  the  ymn  in  (3.3)  may  be  expressed 
in  matrix  form  as  in  Section  2,  of  course.  Let  Am(t0)  be  the  element  signal  vector 
at  time  to, 

-^m(^o)  =  [^ml(^o))  im2(^o)i  •••5  Xm/f(to)]  ,  (3.14) 

Then  ATm(t0)  contains  the  FFT  input  samples  from  element  m  used  in  (3.3).  Also, 
let  Ym  be  a  vector  containing  the  frequency  domain  FFT  samples  from  element  rn. 

Ym  -  [j/ml  t  ym2)  •  •  ■  i  IJmK  ]  •  (3.15) 


Then  Ym  and  Xm  are  related  by 


where,  from  (3.3),  E  is  the  matrix 
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If  A'(to)  is  the  complete  signal  vector  for  the  entire  array, 


f  *i(M 


^2(^0) 


X(t0)  = 


X\i(to) 


and  Y  is  the  vector  containing  all  the  frequency  domain  samples 


T  is  an  invertible  matrix,  of  course.  The  inverse  of  E  in  (3.17)  is  just  the  inverse 


FFT  in  matrix  form, 
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We  have  thus  shown  that  an  array  with  FFT  processing  is  mathematically  equiv¬ 
alent  to  a  tapped  delay-line  array  with  a  linear  invertible  transformation  between 
the  taps  and  the  weights.  In  the  equivalent  tapped  delay-line  array,  the  number  of 
taps  in  each  delay-line  equals  the  number  of  time  samples  used  in  the  FFTs.  and 
the  intertap  delay  Tq  equals  the  sampling  time  T„.  The  theorem  in  Section  2  then 


shows  that  the  array  with  FFT  processing  will  produce  the  same  output  SINR  as 
the  corresponding  array  with  only  the  tapped  delay-lines.  The  FFTs  can  be  inserted 
or  omitted  with  no  change  in  performance. 

An  important  conclusion  that  follows  from  this  result  is  that  FFT  processing 
in  and  of  itself  does  not  offer  any  improvement  in  array  bandwidth  performance. 
The  same  bandwidth  performance  can  be  obtained  simply  by  storing  K  samples  of 
each  element  signal  and  then  weighting  these  samples  directly.  Including  the  FFTs 
between  the  samples  and  the  weights  merely  adds  to  the  computational  burden,  but 
does  nothing  for  the  bandwidth  performance. 


4.  ADDITIONAL  COMMENTS  ON  FFT  PROCESSING 


In  this  section  we  discuss  a  few  additional  points  of  interest  concerning  FFT 
processing. 

4.1  Optimal  Weights  With  and  Without  FFT  Processing 

First,  we  consider  how  the  optimal  weights  with  FFT  processing  compare  to 
those  with  tapped  delay-line  processing.  Let  U  and  W  be  the  optimal  weight  vectors 
with  and  without  the  FFT  transformation  T  in  the  array,  respectively.  Then,  from 
(2.21),  U  and  W  are  related  by, 

U  =  [Tt\~1W.  (4.1) 

However,  the  matrix  T,  given  in  (3.21),  is  symmetrical,  because  E  in  (3.17)  is 
symmetrical.  Hence  (4.1)  simplifies  to 

U  =  T'lW.  (4.2) 

Moreover,  U  and  W  may  each  be  expressed  in  terms  of  element  weight  vectors  Um 
and  Wm  as  in  (2.5).  Eq.(4.2)  is  then  equivalent  to 

Um  =  E~lWm,  1  <  m  <  M.  (4.3) 

Hence  the  optimal  weight  vector  behind  each  element  with  FFTs  is  just  the  inverse 
FFT  of  the  optimal  weight  vector  without  FFTs.  Note  that  because  (4.2)  holds,  the 
same  array  output  signal  is  obtained  with  or  without  the  FFTs, 

YtU  =  \TX\t{T~1W\  =  X  TW.  (4.4) 

as  the  theorem  in  Section  2  requires. 


4.2  Covariance  Matrix  Eigenvalues 


Next,  we  consider  the  eigenvalues  of  the  covariance  matrix  seen  by  the  adaptive 
array  processor  with  and  without  the  FFTs.  These  eigenvalues  are  of  interest  be¬ 
cause  they  control  the  transient  behavior  or  convergence  properties  of  the  algorithm 
used  to  adapt  the  weights  [  1 1] .  Without  the  FFTs,  the  signal  vector  is  X  and  the 
covariance  matrix  is 

<t>x  =  E\X'Xt\.  (4.5) 

Suppose  has  orthonormal  eigenvectors  eXt  and  eigenvalues  A*tl, 

$*6*.  =  \XieXi,  1  <  i  <  KM.  (4.6) 

The  eigenvectors  eXf  satisfy 

=  <5.;,  1  <  i,j  <  KM.  (4.7) 

where  the  superscript  f  denotes  conjugate  transpose  and  <5,,  is  the  Kronecker  delta. 

Now  define  new  vectors 

e„.  =  \/KT~l eZi,  1  <  i  <  KM.  (4.8) 

These  ey<  also  form  an  orthonormal  set.  From  (4.8),  we  have 

elev,=elVKlT-l}'VKT-letj.  (4.9) 

But  because  T  is  symmetrical  (Tt  =  T)  and  71-1  =  ~T‘  (see  (3.23)).  we  have 

\T~1}'  =  =  ±T,  (4.10) 

1  is  a  positive  definite  Herniitian  matrix,  so  it  lias  a  complete  set  of  eigenvectors  and  its 
eigenvalues  are  all  real  and  positive. 


^  -S3 


so  (4.9)  reduces  to 


sleu.  =  KetJ-TT-'e x 


=  ele*, 


=  f>ij,  1  <  i,j  <  KM. 


(4.11) 


Now  in  (4.6)  let  us  substitute  er,  =  ^TeVi  (the  inverse  of  (4.8))  and  multiply 
on  the  left  by  T~1.  We  obtain 


T  l$xTeVi  =  \x,evr 


(4.12) 


Then  we  replace  T~l  by  ^T*  and  T  by  Tr, 


[r**,rK.  =  K\Xievr 


Finally,  from  (2.18)  we  note  that 


(4.13) 


T*xTt  = 


so  (4.13)  is  just 

$ye„,  =  KXx.eVi,  1  <  i  <  /CM. 


(4.14) 


(4.15) 


Eq.(4.15)  shows  that  eu.  and  KXZi  are  the  ith  eigenvector  and  eigenvalue  of  <$>„. 
Thus,  each  eigenvalue  of  is  simply  K  times  the  corresponding  eigenvalue  of  $r. 

From  this  it  follows  that  and  have  the  same  eigenvalue  spread.  (The 
eigenvalue  spread  of  a  matrix  is  the  ratio  of  its  largest  to  its  smallest  eigenvalue.) 
Hence  typical  problems  caused  by  eigenvalue  spread,  such  as  long  convergence  times, 
roundoff  errors  in  covariance  matrix  inversions,  etc.,  will  be  the  same  with  or  without 
the  FFTs. 
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4.3  Weight  D'  anic  Range 
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Now  we  consider  an  issue  of  practical  interest:  how  FFT  processing  affects  the 
dynamic  range  of  the  weights.  We  may  gain  insight  into  this  question  els  follows. 

McClellan  and  Parks  [12]  have  studied  the  eigenstructure  of  the  FFT  transfor¬ 
mation  matrix.  From  their  results,  it  is  eajsily  shown  that  the  matrix  E  in  (3.17) 
has  a  complete  set  of  orthonormal  eigenvector0  e£jt  1  <  j  <  K,  and  that  every 
eigenvalue  A^  of  E  has  one  of  the  four  values  \^K,  —\fK,  +jy/K,  or  —j  s[K.  The 
multiplicity  of  each  eigenvalue  varies  with  K,  the  order  of  E. 

From  the  eigenvectors  and  eigenvalues  of  E  one  can  obtain  the  eigenvectors 
and  eigenvalues  of  T  in  (3.21)  in  an  obvious  way.  Each  eigenvector  of  T  will 
have  (M  -  1  )K  components  equal  to  zero  and  K  components  consisting  of  one 
eigenvector  tEt  of  E.  The  eigenvalues  A^  of  T  will  be  the  same  as  those  of  E  but 
with  multiplicities  M  times  higher. 

Suppose  W  is  the  optimum  weight  vector  without  FFTs  and  U  the  optimum 
weight  vector  with  FFTs.  Then  from  (2.21)  we  have 

U  =  \Tt}-'W  =  T~XW.  (4.16) 

The  optimal  weight  vector  W  may  be  expressed  in  terms  of  its  components  along 
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each  of  the  eigenvectors  of  T, 

MK 

W=Y.aier,-  (4-17) 

;  =  i 

where  each  a;  is  a  scalar  constant.  Moreover,  the  matrix  T~x  may  be  written  in 
terms  of  the  eigenvectors  erJ  and  eigenvalues  A t  of  T  using  the  spectral  decompo¬ 
sition  formula, 
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(4.19) 


Substituting  (4.17)  and  (4.18)  into  (4.16)  then  gives 


MK 


u  =  r^eT. 

h  ^  1 

We  cannot  determine  the  magnitude  of  the  components  of  U  exactly  without  consid¬ 


ering  specific  cases,  because  some  A ?  are  real  and  some  are  imaginary.  But  because 
| At. |  =  \[K  for  every  j,  (4.19)  shows  that  the  weights  with  FFTs  will  generally  be 
smaller  than  the  weights  without  FFTs  by  a  factor  of  about  -if. 


4.4  Performance  Differences  between  Tapped  Delay-Line  and  FFT  Pro¬ 
cessing 

Finally,  let  us  consider  how  the  delay  parameters  are  usually  chosen  in  tapped 
delay-line  and  FFT  processors.  In  the  Introduction  we  noted  that  array  performance 
can  be  poorer  for  FFT  processing  than  for  tapped  delay-line  processing.  Figure  1.6, 
showing  typical  results  with  FFTs,  was  compared  with  Figure  1.4  for  tapped  delay¬ 
lines. 

However,  it  is  clear  from  the  results  of  Sections  2  and  3  that  the  SINR  perfor¬ 
mance  of  an  array  with  FFTs  must  be  identical  to  that  of  the  equivalent  tapped 
delay-line  array.  The  equivalent  tapped  delay-line  array  has  the  same  number  of 
taps  as  the  number  of  samples  in  the  FFTs  and  has  an  intertap  delay  equal  to  the 
FFT  sampling  interval. 

The  performance  difference  noted  in  the  Introduction  is  due  entirely  to  the  fact 
that  typical  comparisons  have  assumed  different  intertap  delays  or  different  numbers 
of  taps  for  the  two  types  of  arrays.  For  a  tapped  delay-line  processor,  the  intertap 


delay  is  often  assumed  to  be  a  quarter  wavelength  at  the  carrier  frequency.  For 
an  FFT  processor,  the  sampling  time  is  usually  chosen  so  the  period  of  the  FFT 
frequency  response  approximates  the  signal  bandwidth.  These  two  amounts  of  delay 
are  usually  very  different. 

Let  us  consider  these  time  delays.  First,  suppose  the  signal  carrier  frequency  is 
w0-  The  time  delay  required  to  produce  a  90°  phase  shift  at  the  carrier  frequency 
(a  quarter-wave  delay)  is  then 


In  general,  suppose  we  have 


(4.20) 


7q  —  rTgo<>, 


(4.21) 


where  r  is  the  number  of  quarter-wave  delays  in  To.  Although  there  is  actually  no 
fundamental  reason  to  do  so,  with  tapped  delay-line  processing  it  is  common  [4j  to 
assume  r  =  l.2 

Now  consider  the  choice  of  T,  in  an  FFT  processor.  We  may  view  the  FFT  in 
(3.3)  as  being  equivalent  to  a  filter  bank.  The  input  to  the  filter  bank  is  xml(t)  and 
the  outputs  are  ymi,  ym 2,...,  ymK •  One  filter  produces  the  output  yml,  another 
produces  the  output  ym2,  and  so  forth.  Let  us  consider  the  transfer  function  of  each 
of  these  filters. 

Suppose  the  input  signal  xmi(f)  is  a  sinusoid  at  frequency  w, 

Xmi(«)  =  e'wt.  (4.22) 

Then,  from  (3.1), 

2 It  is  shown  in  (7|  that  any  choice  of  r  in  the  range  0  <  r  5  ^  will  work  just  as  well. 


(4.23) 


Zmk{t)  - 


eiu,[«-(fc-l)T.] 


For  a  specific  n,  the  output  signal  ymn  may  be  found  by  substituting  (4.23)  into 
(3.3), 

ymn  =  £xmk(t)e-’^k-^ 

k= 1 
K 

—  ^  e->(*-l)lwT,+^-(n-l)le>wt^  (4-24) 

k=\ 

This  may  be  written, 


ymn  =  Hn(  w)e*“*, 


(4.25) 


where  Hn(u),  the  nth  transfer  function,  is 

K 

I 

*=i 

e 


Hn{u)  =  ^e-;(*-»)[^.+^(n-,)| 


_  c-yi^fu)r.+^-(n-i)i'stn  2  [uTt  +  K  (n  1)] 

^  (n  -  1)] 


(4.26) 


In  general,  Hn(ui)  is  a  periodic  function  of  frequency  with  peaks  at  frequencies 

2n , ,  n  —  1,  .  . 

w  =-[»-—].  i  =  -2,  -1,  0,  1,  2,  ...  (4.27) 


For  a  given  n,  the  peaks  of  Hn{u)  occur  every  ~  along  the  frequency  axis.  For 
adjacent  n,  the  peaks  are  separated  by  -j^r.  Figure  4.1  shows  a  typical  set  of 
Hi  (w),  . . . ,  Hk{u)  over  part  of  the  frequency  axis. 

In  studies  of  FFT  processing,  it  is  common  to  choose  T,  so  one  complete  set  of 
K  filter  passbands  approximately  covers  the  signal  bandwidth.  (This  choice  seems 
sensible,  since  it  divides  the  signal  bandwidth  into  K  subbands.)  If  the  signal 
bandwidth  is  Aw,  we  set 


Figure  4.1:  The  Transfer  Functions  #i(w),  Hk{u) 


or 


Ta  = 


27T 

Aw ' 


However,  (4.29)  may  be  rearranged  into  the  form 
T  .  *  w°  4 

I  •  —  **  -  "T -  —  90° , 

2w0  Aw  B 

where  B  is  the  relative  bandwidth  of  the  signals, 

Aw 


B  = 


Wo 


(4.29) 


(4.30) 


(4.31) 


and  Tgo°  is  given  in  (4.20).  If,  for  example,  the  signal  has  a  1%  relative  bandwidth 


Note  that  this  choice  corresponds  to  r  =  400  in  (4.21),  i.e.,  to  an  intertap  delay  of 
400  quarter- waves  in  the  equivalent  tapped  delay-line  array! 

In  Figure  1.4  of  the  Introduction,  which  shows  the  SINR  of  a  2-element  tapped 
delay-line  array,  there  are  two  taps  per  element  and  one  quarter-wave  delay  between 
taps.  In  Figure  1.6,  which  shows  the  SINR  for  a  2-element  array  with  FFT  pro¬ 
cessing,  the  sampling  time  has  been  selected  according  to  (4.29).  Thus,  the  tapped 
delay-line  array  in  Figure  1.4  and  the  FFT  array  in  Figure  1.6  are  not  equivalent. 
The  difference  in  their  performance  is  due  to  the  difference  between  To  and  Tf,  as 
well  as  K. 

An  earlier  report  [7]  discusses  in  detail  how  the  number  of  taps  and  the  amount 
of  intertap  delay  affect  the  SINR  performance  of  a  two-element  tapped  delay-line 
array.  In  particular,  it  was  shown  there  that  setting  Ts  according  to  (4.29)  (i.e., 
setting  r  =  4/B )  makes  r  too  large  to  obtain  optimal  SINR  from  the  array.  For 
optimal  SINR,  r  must  be  in  the  range  0  <  r  <  ~.  Thus,  although  it  seems  intuitively 
sensible  to  choose  T,  so  the  signal  bandwidth  is  divided  into  K  subbands,  in  fact 
this  choice  yields  suboptimal  SINR.  Better  performance  will  be  obtained  if  Ts  is 
chosen  so  the  FFT  period  is  at  least  four  times  the  signal  bandwidth.  The  reader 
is  referred  to  [7]  for  further  discussion  of  this  question. 


5.  CONCLUSIONS 


In  this  report  we  have  shown  that  the  SINR  performance  of  an  adaptive  array 
with  FFT  processing  is  identical  to  that  of  an  adaptive  array  with  equivalent  tapped 
delay-line  processing.  In  the  equivalent  tapped  delay-line  processing,  the  number 
of  taps  in  the  delay-lines  is  equal  to  the  number  of  samples  used  in  the  FFTs,  and 
the  delay  between  taps  is  equal  to  the  delay  between  samples  in  the  FFTs. 

In  Section  2,  we  showed  that  inserting  a  linear  invertible  transformation  between 
the  delay-line  taps  and  the  weights  in  a  tapped  delay-line  array  has  no  effect  on  the 
array  output  signal  or  SINR.  Whatever  changes  are  caused  in  the  signals  by  the 
linear  transformation  are  compensated  for  by  corresponding  changes  in  the  weights. 
Then  in  Section  3  we  showed  that  using  FFTs  behind  each  element  is  mathematically 
equivalent  to  using  tapped  delay-lines  followed  by  a  linear  transformation  of  the 
signals.  From  the  results  of  Sections  2  and  3,  the  main  conclusion  follows. 

In  Section  4,  we  also  discussed  the  effects  of  FFTs  on  the  optimal  weights,  the 
covariance  matrix  eigenvalues,  and  the  dynamic  range  of  the  weights.  Finally,  the 
reasons  for  the  performance  differences  between  FFT  processing  and  tapped  delay¬ 
line  processing  noted  in  the  Introduction  were  discussed. 

The  most  important  conclusion  that  follows  from  these  results  is  that  FFT 
processing  in  and  of  itself  does  not  offer  any  improvement  in  array  bandwidth 
performance.  The  same  bandwidth  performance  will  be  obtained  by  simply  storing 
K  samples  from  each  element  signal  and  then  weighting  and  combining  these  samples 
directly.  Performing  FFT  calculations  between  the  samples  and  the  weights  adds 
nothing. 
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